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Abstract 

A new formalism to calculate electronic states of vacancies in diamond has been 
developed using many-body techniques. This model is based on prevoius molecular 
models but does not use configuration interaction and molecular orbital techniques. 
A generalized Hubbard Hamiltonian which consists of all electron-electron interac- 
tion terms is calculated on the atomic orbital bases. Spatial symmetry T^ and spin 
information of system are included in the form of Hamiltonian, so the eigenstates 
have automatically the correct spin and symmetry properties. Optimizing two key 
parameters of the model that justifies already reported semi-empirical values can pre- 
dict accurate values of the famous absorption lines in neutral and charged vacancies 
i.e. GRl and NDl. With these parameters the location of the low lying ^Ti state 
is 113 mev above the ground state. In addition to these levels whole of the energy 
states of the system is predicted. Since the results are obtained without configuration 
interaction the model can gives the exact contribution of electronic configurations in 
the ground and excited states of the neutral and charged vacancies. The new results 
of the model for the ground and excited states of GRl are reported. 

PACS numbers: 61.72.Bb, 61.72.Ji 



1 



I. INTRODUCTION 



Vacancy is one of the simplest and familiar intrinsic point defects in semi- 
conductors and its effects have attracted high technological and theoretical 
interests. A vacancy in diamond in contrast to the situation in most semicon- 
ductor or insulator crystals is stable in room temperature.^ Since the formation 
energy of vacancies in diamond is high (6 — 7 eV) ^'^ they have been studied 
mainly in irradiated diamond. Optical absorption studies arc used for deter- 
mination of electronic excitation energies of vacancies^'® and EPR, ENDOR 
experiments have been used for measuring the spin and spatial symmetry of 
electronic states. The most famous optical absorption lines in irradiated di- 
amond are GRl (with zero phonon line at 1.673 eV)^^ ^nd NDl (with zero 
phonon hne at 3.149 eV)^^. The former is attributed to a neutral vacancy-*^^ and 
the latter to a negatively charged one.^^ 

In addition to these two famous lines many other absorption lines have been 
reported such as: N3 (2.985 eV)^ H3 (2.463 eV)^^ NV (1.945 eV)^^ which 
have been attributed to vacancy- nitrogen^^'^°, also GR2 to GR8 hnes^ which 
have been attributed to the natural-charged vacancy complexes.^^ Despite the 
simplicity of the system, different models have been suggested to explain the 
properties of these optical absorption lines. Generally there have been two 
main approaches to this problem, The first group of approaches are localized 
models^^"^^ and the second ones are extended models. In the first group, 
the electrons of the broken bonds which are in the tetrahedrally symmetric local 
potential of the vacant site are considered as an isolated molecule (many body 
approaches). In the second group i.e. extended models, the effect of lattice 
is more important on the vacancy and it is considered as a missing atom in a 
supercell with few hundreds of atoms (single particle approaches). Similarity 
between these two categories is their group theoretical aspects for accounting 
symmetries of the system while their main difference is in the importance of e-e 
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correlation effects. Lannoo and Bourgoin ^^liave developed a model with four 
parameters which briliantly explaines the validity range of the single and many 
particle approaches for this problem. According to Messmer and Watkins^^ a 
successful model for describing vacancy energy states should explain the fol- 
lowing points: The position of the defect energy levels relative to the valence 
band edge, the electronic wavefunction of defect electrons for comparing with 
experimental data, lattice relaxation and distortion in neighborhood of vacant 
site and it must be the basis of a practical computational scheme. The localized 
approaches have relatively successful answers to these questions. 
The first suggested model in this regard is the famous model of the Coulson 
and Kearsly^^. Their model suggests a computational scheme that considers 
e-e interaction in the system accurately. They used group theoretical methods 
to manually construct symmetry and spin adapted wavefunctions for including 
spin and symmetry considerations of vacancies. By applying configuration in- 
teraction technique to the calculation results, they predicted qualitatively good 
results for electronic states of the vacancies, their model used two semi-empirical 
parameters to obtain satisfactory result for energy of the GRl transition. By 
using these semi-empirical parameters their predictions for NDl transiton and 
also the posision of the low lying ^Ti state disagrees with avaiable experimen- 
tal data. Up to now this model has been the most successful computational 
scheme for describing physical properties of the optical absorption lines of va- 
cancies in diamond, ^'^^ and the next suggested localized models^^"^^ are based 
on the framework of this model. 

Localized models, unlike the nowadays commonly used density functional the- 
ory (DFT) methods which pay more attention on the ground state information, 
can also give the excitation energies of the vacancies beyond the one electron 
approximation. 

For improving the quantitative agreement between theory and experiment we 
introduce a new approach for localized models. This approach is based on a gen- 
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eralized Hubbard Hamiltonian for electrons of the vacancy system with atomic 
orbital bases and uses computer adaptable many body techniques instead of 
usually used molecular orbital method (fourth item of messmer and watkins). 
In this paper at first we review the computational scheme of the present model 
and discuss the advantages of the new notation which is used for this problem. 
Then we apply this new scheme to solve generalized form of Hubbard Hamil- 
tonian for neutral and charged vacancy systems. For evaluation of solutions of 
this Hamiltonian we start with parameters which were already reported. Our 
using value for the two semi-empirical parameters slightly differs from already 
rpoeted values but other six parameters are as the same as the theory calcula- 
tion. We report the effect of the hopping parameter on the energy spectrum of 
this Hamiltonian by varying this parameter. 

The model predictions for GRl and NDl transitions with justified parameters 
and comparison the results with experimental data will be discussed. The lo- 
cation of the low lying ^Ti state is predicted and we disscuss its comparision 
with other theoretical reports and experimental evidences. The results of our 
model can justify the semicmpirical values that already were reported. '^'^'^-'^ We 
also will show that, this model is independent of familiar configuration interac- 
tion technique and the electronic configurations mixing of each state comes out 
in a natural manner from the formalism. The model predicts whole of the ex- 
act eigenstates and enery levels of the systems with the justified semi-empirical 
values. Finally the new information concerning electronic configurations of the 
ground and excited states of GRl transition are reported which can be used by 
other approaches to this problem. 

In our calculation the effects of lattice relaxation were ignored, which according 
to recent theoretical works are small.^^"^^ 
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II. COMPUTATIONAL MODEL 



The present computational model takes into account electron correlation ef- 
fects accurately and gives the correct spin and spatial symmetry of eigenstates. 
It is also suitable for small atomic cluster systems or molecules provided the 
hopping, Coulombic and exchange integrals are known. In the localized model 
framework, it is assumed that physical properties of neutral (charged) vacancy 
depend on four (five) electrons of adjacent dangling bounds. The exact Hamil- 
tonian of such a system can be written as Eq. (1). A more simplified form of 
this Hamiltonian recently is applied to carbon nanotube systems ^^'^^ . 

iif ijlmaa 

Where i,j,l,m. are atomic sites indices which go from 1 to 4 , cr is spin index 
which is —1/2 or +1/2 , cj^ and Cia are operators which create and annihilate 
electrons with specified spin a on site i. tij and Vijim parameters are hopping 
and e — e interaction overlap integrals which are defined as follows: 




in the above integrals ijji is single electronic wave function in a specified site. 
The advantage of using second quantized form of Hamiltonian for this prob- 
lem is its capability to include the spatial symmetry and spin information of 
the system in the form of Hamiltonian, hence the eigenstates have automati- 
cally the desired spin and symmetry properties. Also the effect of the changing 
each parameter directly can be observed on the electronic energy levels and 
also wave functions. This new approach uses atomic orbital bases to solve the 
Hamiltonian. In the previous localized models it was common as a start point, 
to manually construct symmetry and spin adapted molecular wave functions as 
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correct basis for accounting symmetry considerations of the system then these 
basis should be used in evaluating Hamiltonian matrix elements. This process 
is a laborious task and can not be incorporated in a computational software 
easily. 

Unlike previous methods, present computational scheme needs significantly less 
effort and is straightforward in converting to computer software language. The 
starting point of the model is constructing appropriate many-body basis needed 
for solving the Hamiltonian. 

All of the possible distribution of electrons in the four adjacent atomic orbital 
make a set of complete basis in the configuration space. In contrast to previous 
models that have used molecular orbital for such a system we have developed 
atomic orbital bases by taking into account z component of the spin degree of 
freedom with occupation condition of each orbital. Following notation can be 
used to represent each configuration state. Ignoring the spin degree of freedom 
simply reduces the states to the occupation number or the Fock space repre- 
sentation. Fock space representation is commonly used in modern molecular 
calculations.'^^ 

I^i >= |aio,aa,ai2, •••,ai8 > (4) 

Parameters Oji to ajg are -1 or -|-1 (for spin down and up on the site i respec- 
tively) which show occupation condition of each state. In fact these coefficients 
not only show extension of electronic wave function on each neighboring site of 
vacancy but also include spin information of the system, a^o is the sum of all 
aj's i.e. the total spin of each along z axis. We will show the importance 
of this quantum number in our computation scheme later. Other capability of 
this notation is estimation of the maximum dimension of configuration space 
which one needs for calculating the Hamiltonian for In the 

system there are four (five) electrons respectively in four orbital, so the possible 
configurations of these electrons in the 8 accessible states (2 spin and 4 space 
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TABLE I: In this table, Sf is total spin in z axis direction, Hup is number of electrons 
with spin up and n-q, is number of states with specific Sz- 

or site degree of freedom) according to statistical combination rules are: 

Ct = 70 (V) (5) 

Cl = 56 {V-) (6) 

The fractions of 70 (56) configuration states in {V~) system which belong to 
each Sz block are summarized in Table I. This counting method is in contrast 
to the old molecular orbital approaches which one needs to allow construction 
of symmetry and spin adapted basis by hand and then counting them. 
After constructing the set of complete basis, we will calculate two form of Hamil- 
tonian in this space. The first one is approximate form of Eq. (1), i.e. a simple 
extension of Hubbard model and the second is the exact form of the Hamiltonian 
(Eq. (1)) which we call it generalized Hubbard Hamiltonian. 

III. GENERALIZED HUBBARD HAMILTONIAN 

The main goal of the Hubbard Hamiltonian^^ is resuming atomistic nature 
of the solid besides the free electron gas theory. This model assumes that the 
most important part of e — e interaction terms is on site term on the Columbic 
parts of Hamiltonian. The conventional form of this Hamiltonian consists of 
two parts, hopping term or t term and on site term or U term which is Vun 
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term of Eq. (3). 

-H" = tijcl^Cja + UiUi t rij i (7) 

ija i 

Famous Hubbard model with two parameters t and U is not sufficient for solving 
the vacancy problem, the results of our model with these two parameters are 
not satisfactory. This result is physically expecting since assuming only two 
parameters of the Hubbard model for this system means that the electrons 
interact with each other only when they are at the same site of at the same 
atom of the vacancy and this eliminates the interaction of the electrons when 
they are in the different atomic site of the vacancy. Therefore we started with 
extended Hubbard Hamiltonian with one extra parameter V which is Vijij term 
of Eq. (3). This Hamiltonian is equivalent to the Lannoo modeP^ where the U 
parameter of their model is called V in our model. This means that the on site 
and two site direct overlap integrals are the most dominant ones and the rest 
parameters are negligible (Eq.(8)). 

H tijc\^Cja + ^ UiUi ]nii+]^ VijUiaUj^' (8) 

ija i i^jaa' 

Since the atomic sites in vacancy are from the vertex of a symmetrical tetra- 
hedral, the distance between each two of them is the same, and hopping (ij^), 
on site {Ui) and two site iVij) parameters of Hamiltonian are all independent of 
i,j indices and can be put outside of the sum Eq. (8). 

H c\^Cj^ + U^ni] Hi i niaUj^' (9) 

As discussed in sec. H, for system the size of the set which includes 

complete basis for solving this Hamiltonian is 70 (56). As a result solving Hamil- 
tonian in this set results in a 70 x 70 (56 x 56) Hamiltonian matrix for 1/° (V^) 
respectively. As the starting point we used numerical values of Hamiltonian 
parameters t, U and V which were reported by semi-empirical and theoretical 
works on diamond vacancies. ^^'^^ 
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This formalism can investigate the effect of variation of each parameter on the 
energy spectrum directly. By this means we found that variation of parameter 
t has not any effect on the GRl transition energy. But variation of parameters 
U and V changes the transition energy of GRl. These variations do not affect 
the sequence of the levels. The results of calculation for show that for a 
wide variation range of variables, t greater than —12.5 eV, U greater than 8.2 
eV and V greater than 4.0 eV, the ground state of the system has symmetry 
(double degenerate spin-less) and we can obtain appropriate transition to an 
excited state (triple degenerate spin- less). With our model we were able to 
obtain the experimental value of 1.673 eV for the GRl transition. 
For system, variation of t has a dramatic effect on NDl transition energy, 
but does not alter the ground and excited states sequence. Unlike this, U and V 
parameters which have only small effect on the values of the transition energies 
of the system. For the same range of variation of parameters we had a transition 
from '^A2 ground state to a ^Ti excited state {NDl transition) but By the same 
parameters which we obtained the transition energy of GRl (1.673 eV), there 
is a wide gap between the ground and excited states of NDl transition (about 
30 eV). 

Only with parameter t close to —1 was the result close to the experimental 
value. This high value for t seems unreasonable for the vacancy problem and is 
very far from reported range.^^'^^ 

These results of the extended Hubbard Hamiltonian are in agreement with the 
results of the Lannoo and Bourgoin model which has simplified the full Hamil- 
tonian in the first quantization form by four parameters. They have used sym- 
metrical molecular orbital base and also configuration interaction similar to 
Coulson and Kearsly model. The benefith of the reduction of the number of 
Hamiltonian parameter in their model was in this point that the validity range 
of single and many particle approches could be investigated. They investigated 
the results of their model in a full range of parameters and also used the semi- 
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TABLE II: List of estimated Hamiltonian parameters in Hubbard model and diamond 
vacancies problem. 

empirical parameters of the Coulson and Kearsly model. They successfully 
explained where the many body model works and where the single particle ap- 
proaches can be valid. Lowther^^ applied a similar model to the Lannoo model 
to explain some related optical phenomenon of the vacancies in diamond. Main- 
wood and Stonham^^ applied Lannoo model directly to the vacancies electronic 
states in diamond. It seems that lack of the quantitative agreement with ob- 
served NDl transition energy (3.15 eV) from ground state to a ^Ti excited 
state in V~ system is due to elimination of exchange terms in Hamiltonian of 
vacancy system. Since calculation results of V~ system show weaker agree- 
ment with the experiment than those for V^, we can deduce that the role of 
the exchange terms in charged electron vacancy which has higher electron den- 
sity is more important than that the four electron (neutral) one. Quantitative 
failure of the results of this simple extension of Hubbard Hamiltonian reveals 
that for vacancy system e — e correlation effects are important and can not be 
neglected in contrast to the other cases which we are able to apply the Hub- 
bard model. This point also was emphasized by previous models which have 
attempted diamond vacancy problem. ^^^^^''^^''^^ Comparing numerical values of 
overlap integrals of Hubbard model^° and diamond vacancy problem^^'^^ (Table 
II) we find that the values of exchange integrals in this problem are comparable 
to the values of direct integrals so their elimination are not reasonable. 
This leads us to keep all of the, t, U and V terms besides exchange terms in 
four (neutral) and five (charged) electronic Hamiltonian of the system similar to 
Coulson and Kearsly model. ^^'"^^ We will call it generalized Hubbard Hamilto- 
nian (Eq. (1)). Although this significantly increases the volume and complexity 
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of the calculation, but since we take into account exact e — e correlation effects, 
it makes our model more realistic. 

We return to Hamiltonian form of Eq. (1). Similar symmetrical argument is 
valid for a system with symmetry and we can put hopping parameter {t) out 
of the sum (Eq. (10)). 

H ^t^c\^Cj^ + ]^ Vijimclc^^^Cma'Cla (10) 

Symmetry of defect molecule, reduces the number of the independent parame- 
ters Vijimi from 4^ (256) components to only 7 independent ones, hence symme- 
try considerations can be included in the generahzed form of Hamiltonian with 
a total of 8 parameters. These consist of the hopping t, two direct Coulombic 
integrals (C/, V) and five exchange {XI, ...,Xb) integrals as follows: 
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Since the Hamiltonian is spin independent, useful conserving quantum numbers 

[H,S']^[H,S,]^0 (11) 

arc total spin of four electron system (S^) and the z component of total spin 
(Sz). The important computational point about Eq. (11) is that it should hold 
for each arbitrary set of parameters of Eq. (10). 
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By these conservation rules we can easily decrease dimension of the nonzero 
block of Hamiltonian matrix by using conservation of total spin (S^) and z 
component of the total spin (3^). Therefore the volume of calculation signifi- 
cantly reduces and matrix element evaluation only confines in subspaces where 
wave functions have the same values of Sz- As a result, the maximum dimension 
of nonzero matrix block in the case of (V^) systems becomes 36 (24) respec- 
tively. The final step in solving the Hamiltonian of the (^~) for obtaining 
the resultant eigenstates with definite Sz and S'^ values is the transformation of 
Hamiltonian matrix from to {S'^, S^) basis. 

We have used the parameters which were calculated by Coulson and Larkins^^ 

to investigate energy spectrum and eigenstates of the system. These parame- 
ters have been calculated using Slater and Clementi type functions as atomic 
orbitals. 

The parameter calculation is independent from the Hamiltonian solution formal- 
ism. Since the parameters are calculated two times by the Coulson ^^'^^ during a 
long period, other models similar to us have not attempted to re-estimate their 
values and all of them have used the results of the Coulson estimation for the 
parameters. 2i,22,26,32r]n]-^g calculation of parameters were based on symmetric 
and antisymmetric molecular orbital basis (a, tx, ty, t^jwhich were combinations 
of a, 6, c, d atomic orbitals. These atomic orbitals belong to the four nearest 
neighbor atoms of vacant site and are oriented inward to its center, however 
in the generalized Hubbrad Hamiltonian Eq. (1), parameters are calculated di- 
rectly by using overlap integrals of these a, b, c, d atomic orbitals. These overlap 
integrals can be extracted from the older symmetric and anti-symmetric ones 
by changing the calculation basis. Related calculations have been performed 
and the resultant parameters are used in the computational scheme. 
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TABLE III: List of eight parameters that are needed to evaluate energy states of 
vacancy. They are estimated by this work and other different calculations. 

IV. RESULTS AND DISCUSSION 

The two key parameters of this model are t and U . These two have been 
evaluated semi-empirically in the previous models and were starting point in 
calculation of the other parameters. ^^'^^ The reported values for hopping pa- 
rameter t cover a wide range of variations, for example both —7.13 eV and 
— 16.34 eV values at the same time are reported for this parameter. ^^'^^ Con- 
cerning parameter U , the semi-empirical value is estimated to be 13.29 eV and 
after revising 12.85 eV while the theoretically calculated value is 19 eV .^^'^^ 



In evaluation of the t parameter despite of other 7 Coulombic parameters, one 
not only need the information of the snatial distribution of the wave functions 
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FIG. 1: Variation of ground and excited states of GRl transition, with respect 
to hopping parameter t. The t = —6.25 eV, —10.22 eV are changing point of the 
ground state with ^Ti state. 
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but also the form of the psedupotential of carbon atoms is needed, so the uncer- 
tainty in the value of this is more than other parameters. In the present model 
the variation of the energy spectrum with t has been investigated and results 
are shown in Fig. 1. The values of the other parameters are according to table 
III. According to this figure, for obtaining correct results for ground state of 
i.e. ^E, parameter t should be confined in the following range. 

-6.25eV > t > -10.22e1/ (12) 

Choosing t out of this range converts the sequence of two lowest states i.e. 
and ^Ti, and gives wrong ground state ('^Ti) for the system. Increasing value 
of the t will lower the experimentally observed spin quintet state. It is 
interesting to note that above the upper limit of t, ^A^ competes with ^Ti to 
be the ground state and at t value equals to —5.96 eV, becomes ground 
state of system. It seems that by increasing t we reach an area where the 
Hund's rule can be applied (according to this rule ^^2 must be the ground 
state). This might be due to the point that t is the indicator of the looseness 
of the electrons from the nuclei so by increasing t we can consider the electrons 
as more delocalized electrons. 

For charged vacancy case, the calculation results are shown in Fig. 2. From this 

figure we conclude that for obtaining '^A2 as ground state for the system, pa- 
rameter t should be restricted as: 

-6.95el/ > t > -WMeV (13) 

The lower limit is the negative of the ionization energy of electron in C-C bound 
of diamond. From this figure we can find that increasing t over upper limit 

converts the ground state to ^Ti which disagrees with the experiment. 

By comparing the results of the calculation for both V'^ and V~ we can put new 
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FIG. 2: Variation of ground and excited ^Ti states of NDl with respect to 
hopping parameters t. The i = —6.95 eV is changing point of the ground state with 
the excited state. 

boundaries to the single particle hopping parameter t as follows: 

-6.95ey > t > -W.22eV (14) 

This restriction on t parameter, which is the most uncertain parameter in the 
vacancy electronic states problem, is tighter than already reported range ( —7.13 
eV to —16.14 eV).^^'^^ Variation of t in this allowed range does not have signifi- 
cant effect on GRl transition energy. This is in agreement with other theoretical 
works. ^^'^^ However this variation has a dramatic effect on NDl transition en- 
ergy. 

Numerical optimization of parameter t in this range and also parameter U gives 
—7.747 eV and 13.58 eV for t and U respectively. These values are very close 
to the semi-empirically derived values by Coulson and Larkins (the difference 
between our optimized value by semi-empirical values in this two parameters are 
less than 8 and 5 percent respectively) . Other six parameters are as the same 
as the theoretically calculated values by them using Slater type functions 
Calculation methods of these parameters are completely independent from the 
Hamiltonian solution formalism so the previous models 21,22,26,32 gjj^iig^r 
have used the parameters which are reported semi-empirically and theoretically 
by Coulson and Kearsly ^^and Coulson and Larkins.^^ These parameters are 
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calculated by them during a 17 year period and with the Slater type functions 
so it seems that the parameters are well justified. The difference between our 
suggested parameters with two semi-empirical parameters of their model is very 
less than the wide range of variation of these parameters in their model. ( pa- 
rameter t from —7.13 eV to —16.34 eV and parameter U which smi-empirically 
is 12.855 eV and theoretically is 19 eV) 

With these justified values the present model can simultaneously obtain the 
exact transition energies of the GRl (1.67 eV) and NDl (3.15 eV) for neutral 
and charged vacancy system. 

The set of our parameters and parameters that were previously computed based 
on Slater and Clementi type functions^^'^^'^^ are listed in Table (III). The model 
predicts the wrong ground state with parameters which are based on Clementi 
type functions, i.e. ^Ti instead of the experimentally observed state for 
This is similar to the Coulson and Larkins's report. However parameters 
which are based on the Slater type functions predict correct spin and symmetry 
for ground and excited states of GRl and NDl absorption lines but the values 
of the transition energies are not satisfactory. The results of the present model 
calculation are summarized in Fig. 3,4. 

For neutral vacancy as it is shown in Fig. 3, there is a transition from double 
degenerate spinless ground state to a triple 

degenerate spinless excited state {^E — >^ T2) with a transition energy equal 
to 1.67 eV. 

We also see that there is another level ^Ti, very close to ^E ground state 
with a very httle energy difference i.e. about 113 meV above it. This value 
agrees with the recent experimental expectations which predict it will be more 
than 100 meV .^^ However the results of Coulson models gives 40 meV which 
disagree with the experiment. "^^'"^^ This very close energy level to ground state 
also has been used to explain the behavior of R2 center absorption line in 
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FIG. 3: Neutral vacancy energy levels with, t = —7.747 eV, U = 13.58 eV (present 
parameters) gives 1.673 eV for to transition. The right-hand side levels are 
the continuation of the left-hand side levels. 

diamond where it was attributed to the strongly perturbed vacancy. By this 
assumption the energy difference of this level with the ground state has been 
predicted to be between 40 meV and 200 meV.^^ Mainwood and Stonham 
have used the simple model of Lannoo and successfully explain many 
physical features of the vacancy. Lannoo for simplifying the problem do not 
enters the exchange parameters into energy levels and use only two parameters 
for describing energy levels. This is unlike to Coulson and our model 
which reserve whole of the e-e terms. Besides explaining many physical feature 
of the vacancy they gave some estimation for the position of this level. Using 
the model of Lannoo they positioned ^Ti at 200 meV above the ^E, Although 
this value is as the same magnitude as the error in the positions of the levels 
in their calculations.^^ 

Our reported location for ^Ti arises in a natural manner similar to GRl and 
A^Z^l after evaluation of the exact Hamiltonian of the system with the justified 
parameters. The error which can be attributed to this value in our calculation 
is the difference between our parameters and semi-empirical parameters of 
the Coulson. This error is less than 8 percent which is significantly less than 
already reported error for the position of this level in similar calculations.^^ 
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In obtaining this value we have not assumed that the vacancy is highly 
perturbed as it already was reported.^^ This value can help EPR experimental 
investigations that expect to observe this state in temperatures higher than 
300 



T 100 



33 - 82.512 'T -98 

82 - 96 



82 - - 94 ^ 

> I 

» 81 - 92 ^ 

— UJ 

LU 2-,- 

91 , 90.286 T „. 

80.350 



80.253 E 

79 2r 87.161 

/3./b3 T 86.520 2j 

79.364 "A 86.272 



FIG. 4: Negative vacancy energy levels with, t = —7.747 eV, U = 13.58 eV (present 
parameters) gives 3.159 eV for to ^Ti transition. The right-hand side levels are 
the continuation of the left-hand side levels. 



In the charged vacancy system, as it is shown in Fig. 4, the model predicts 
a transition from a non-degenerate ground state to a triple degenerate excited 
state with a spin equal to 3/2 (^^2 — >^ Ti). The energy difference is 3.15 eV 
which is the famous observed NDl absorption line. Previous models that eval- 
uate e-e interaction exactly can not obtain this value by any set of parameters. 

33,41 

In summary as it has been mentioned by other authors ^^'^^ the previous models 
that evaluate e-e interaction exactly ^^'^^ can only explain GRl transition en- 
ergy quantitatively. They could not predict NDl transition energy and also 
they have predicted the experimentally wrong energy for the low-lying ^Ti 
level. These unsatisfactory quantitative results are obtained not only by the 
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semi-empirical parameters but also in a wide variation range of these two semi- 
empirical parameters. ^^'^^ This variation ranges is very wider than the difference 
of our optimized values with semi-empirically reported ones for t and U. Other 
six Coulombic parameters in our model and the previous ones are evaluated by 
theory. ^^'^^ Our model can simultaneously give accurate transition energies for 
both of the GRl and NDl transitions besides the experimentally acceptable 
energy for the low-lying ^Ti level with almost the same parameters ( only two 
semi-empirical values differs less than 8 percent). 

The results of our model have justified the Coulson and Kearsly semi-empirical 
parameters and also their guess^^ which expected that the t — —S eV value be 
more realistic than t = —16.34 eV . Both values of this parameter have been 
used in their models. ^^'^""^ 

All of these qualitative and quantitatively satisfactory results in our model have 
obtained without configuration interaction. The good quantitative results for 
GRl and even other reasonable quahtative results in the Coulson and Kearsly 
model and all of the previous molecular approaches 22-26,32 ^^^y be achieved 
with assistance of the Configuration interaction technique. This is for the first 
time that a molecular approach put aside the configuration interaction tech- 
nique and could predict more experimental data. Configuration interaction has 
been widely used in all of the previous theoretical works to overcome the dis- 
agreement between theoretical results and experiment even at the qualitative 
levels. Although precise review of the previous exact models ^^'^^ re- 
sults on the charged vacancy with parameter t = —16.34 eV shows that applying 
configuration interaction gives wrong ground state for V~ . 
Similar ambiguity also exists in the electronic configuration of the ground and 
excited state of the GRl. While the models assume that the ground and excited 
states of GRl arise from the single configuration aH'^, at the same time they 
try to find or use some semi-empirically coefficients for the contribution of the 
other possible configurations in the ground and excited states of Ci?!.^^'^^ Here 
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we refer to Davies ^ remark in this regard: 

"The importance of configuration interaction at the vacancy in diamond was 
controversial, and alternative approaches were investigated in which emphasize 
was on the benefit of using large cluster" . 

The main advantages of the presented formalism which is based on molecular 
model is removing the configuration interaction besides it's satisfactory quanti- 
tative results. By this means our model have solved one important shortcoming 
of the molecular approaches with respect to cluster models or density functional 
approaches. 

In the present model, the resultant eigenstates of the Hamiltonian have unique 
expansion on the starting constructed basis in configuration space, so config- 
uration interaction comes out in a natural manner from the formalism. From 
this point, it is possible to find the contribution of each electronic configuration 
in ground and excited states of the system. 

For GRl transition, the results of calculation for the ground state and 
the excited state with the justified parameters are illustrated in Fig. 5. 

The possible electronic configuration for states with value of S* = are: 
(1, 1, 1, 1), (2, 1, 1, 0), (2, 2, 0, 0) where the numbers in parenthesis show the oc- 
cupation number of each tetrahedral orbital of vacancy system. As it is shown 
in Fig. 5, 

all possible electronic configurations are present in the ground state ^E. The 
most probable configuration in ground state arise from (1, 1, 1, 1) configuration 
(61 percent) which has only singly occupied orbitals. Two doubly occupied 
configurations have very low contribution (3 percent) in the ground state. In 
(1, 1, 1, 1) configuration the electrons have maximum separation distance, which 
results in lowering the Coulombic repulsion energy of the system. In contrast to 
this, the excited state of GRl, i.e. is orthogonal to (1, 1, 1, 1) configuration 
and is almost a purely doubly occupied configuration state. This configuration 
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FIG. 5: Percentage of the contribution of each possible electronic configuration in 
ground and excited state of GRl transition for neutral vacancy system, the numbers 
in parenthesis are occupation number of each tetrahedral orbital 

has the highest Coulombic repulsion energy due to its lowest separation distance 
between electrons of the system. 

Roughly speaking, we can state that GRl transition is resultant of a transition 
from higher probability configuration (1, 1, 1, 1) and (2, 1, 1, 0) configurations of 
the system to (2, 1, 1, 0) configuration . 

Our results modifies the common belief about ground and excited 
states of the GRl transition that assumes they both belong to a^t^ 



As the results of our calculation show, the absence of (1, 1, 1, 1) configuration in 
the ground state of the system is unexpected, since it has minimum Coulombic 
repulsion energy among the other configurations of the system. Also it is the 
only compatible configuration of the unique excited state, which in the de- 
localization limit (t > —5.96 eV) is the ground state of the system. 
For NDl transition, since the ground and exited states have the spin equal to 
3/2 the only configuration which can be encountered is (2,1,1,1). Therefore 
the ground and excited states are only a combination of these configurations 
with different Sz values and the results are same as other models . 
The new information about the ground and excited state of the GRl is 



configuration. 
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very important, It can open a new opportunity for ab-initio density func- 
tional approaches which have not been attempted to estimate the value of this 
transition. This is due to the results of the previous molecular orbital models 
which have assumed the ground and excited state belong to the same configu- 
ration o?t^. However by the new result of present model, they can try to obtain 
the transition energy of the GRl transition similar to NDl. These information 
are also useful for theoretical approaches which start from a guess wave function 
for ground and excited states of such a molecular system. 

V. CONCLUSION 

Theoretical studies of the vacancy electronic states in diamond have been a 
challenging problem for the five decades. Localized models which originate from 
Coulson and Kearsly pioneering work, predict correctly the spatial symmetry 
and spin of electronic states but can not give quantitatively good results for 
electronic transitions in diamond vacancies. In their model despite of excellent 
theoretical framework and strong physical intuition, the computational method 
are very old and can not be adapted with nowadays powerful computer assisted 
methods. 

We have developed a new approach based on generalized Hubbard Hamiltonian 
and many body techniques. A complete set of atomic orbital basis for describing 
this system is constructed according to Sz representation of each electronic 
configuration in the four vacancy orbital. This confines the calculation only in 
definite subspace. This representation can simply estimate number of states 
of the system and also each definite Sz block size of Hamiltonian. This is in 
contrast to previous models that uses molecular orbital bases and manually 
construct each symmetry and spin adapted wave functions. 
The generalized Hubbard Hamiltonian which contains whole of the e — e 
correlation terms is included with spatial symmetry information of the vacancy 



22 



therefore the number of the independent parameters of the Hamiltonian is 
reduced to 8. By this means the effect of the parameters on the energy states 
directly can be examined. A detailed computer software is developed to solve 
this Hamiltonian . The resultant Hamiltonian matrix is transformed to change 
basis from Sz to S"^ representation. As the result, the eigenstates have definite 
and correct spin and symmetry properties. 

The effect of changing each Hamiltonian parameters on the electronic states of 
the system directly can be examined. This calculation shows that for obtaining 
experimentally correct results for the ground state of y° and V~ , parameter 
t should be in the —10.22 eV to —6.95 eV range, although increasing t in 
this range has low effect on GRl but significantly increases NDl transition 
energy. According to this model the optimum value for two key scmi-cmpirical 
parameters t (in allowed range) and U are —7.747 eV, 13.58 eV respectively. 
Two famous absorption lines of and V' namely GRl (1.67 eV) and NDl 
(3.15 eV) can be obtained simultaneously and exactly with these parameters 
in addition to the location of the low- lying ^Ti state at 113 meV above the 
ground state of the V^. These results have been obtained for the time by 
such justifiable parameters. These parameters also justify already reported 
semi-empirical values that only could estimate GRl transition and gave 
experimentally wrong results for ^Ti. 

For the first time the configuration interaction procedure, has not used in such 
a calculations and it arises in a natural from from the formalism. Despite of 
its ambiguous physical nature, the configuration interaction techniques plays 
a fundamental and vital role in reaching to an agreement with experimental 
data in previous localized models. This vital role even exist at the level 
of quahtatively acceptable results. This point enables us to estimate the 
contribution of each electronic configuration in ground and excited states of the 
system. The results for GRl absorption line show that it is almost a transition 
from configurations that consist of a pure singly occupied orbital and a singly 
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double occupied orbital to a singly occupied orbital configuration. This new 
information can be used by the density functional approache to this problem. 
The results of our model justifies the first suggested molecular model by the 
Coulson and Kearsly is the five decades ago and shows that their model has 
enough physical depth to accompany with a new computational methods to 
obtain very satisfactory quantitative results. The results also show the benefit 
of the localized models with respect to cluster or density functional approaches 
by no need of configuration interaction . 
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